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We propose that in noncentrosymmetric superconductors with weakly asymmetric spin-orbit in- 
teraction the field-induced pair correlation between the spin-orbit split different bands ignored in 
previous studies yields unique effects; i.e. the Pauli depairing effect is anisotropic in the momentum 
space, and as a result, magnetic fields induce point-node-like anisotropic gap structure of the quasi- 
particle energy even for isotropic s-wave states, which seriously affects thermodynamic quantities 
at low temperatures. Also, it is shown that when the magnitude of the spin-orbit interaction is 
smaller than the superconducting gap, the specific heat as a function of a magnetic field exhibit 
a two-gap-like behavior, even when there is only a single gap. These features characterize parity 
violated Cooper pairs in weakly noncentrosymmetric systems. We suggest the possible detection of 
these effects in the superconductor with weakly broken inversion symmetry Y2C3. 



I. INTRODUCTION 



One of the intriguing features of the recently discovered noncentrosymmetric superconductorai*2ii2ii^i^ is the realiza- 
r tion of parity violated Cooper pairs, which leads various exotic phenomena.-'^i^i^ii^'iiii^'i^ii^ii^'i^iii'i^ii^i^S'^ii^^ The 

^ ^ investigations for unique effects characterizing parity violation have been mainly focused on the case with strongly 
^ broken inversion symmetry (IS), where the spin-orbit (SO) band splitting ^Ego is enormous compared to the supercon- 
C/3 I ducting (SC) gap A, because it is commonly believed, and partly true, that effects due to broken IS are more prominent 
for stronger SO interaction. However, in the present paper, we propose the possibility of novel phenomena which 
characterize parity violated Cooper pairs in a unique way inherent in weakly noncentrosymmetric systems. These 
phenomena are essentially raised by the Zeeman effect on pairing states in the SO split two bands. In the weakly non- 
centrosymmetric case Ea,o ~ A, the Zeeman magnetic field induces substantial pair correlation between the SO split 
fH different bands competing with the intra-band pairs, which is, in contrast, negligible in strongly noncentrosymmetric 
Q systems with E'so ^ ^^ ,14, 15, 16,18^19,20 demonstrate that the field- induced inter-band pair correlation changes 
O . drastically low-energy properties of the SC state in the case with weakly broken IS, yielding the following unique 
effects; (i) the Pauli depairing effect is anisotropic in the momentum space, and as a result, magnetic fields induce 
the point-node-like structure of the excitation gap even for isotropic s-wave states, which seriously affects thermody- 
^ namic quantities, yielding distinct behaviors of them at low temperatures, (ii) For i?so < A, the specific heat as a 
function of magnetic fields exhibits a two-gap like behavior even when there is only a single SC gap. These effects 
0^ are associated with the momentum-dependent spin orientation of Cooper pairs which characterizes parity violation. 
Thus, our results suggest a possible new direction of the experimental search for parity-violated Cooper pairs. 
. As a matter of fact, our findings are relevant to the recent experimental studies on the weakly noncentrosymmetric 
C ■ superconductor ¥2031^1^1^1^ This system has a cubic crystal structure with the space group symmetry /43c? breaking 
, inversion symmetry. Thus, the asymmetric SO interaction can be approximated by the Dresselhaus type interaction. 
■ The unique feature of this system is that the SO splitting is almost of the same order as the superconducting gap)^i2£ 
r\ \ and thus the situation considered in the present paper may be realized in this material. In the last part of this paper, 
J> ■ we shall compare our results with the experimental observations for this system, and discuss the possible realization 
k>( ' of the distinct phenomena characterizing parity violated Cooper pairs in Y2C3. 

. Although the existence of the indispensable inter-band pairing correlation in addition to the intra-band pairs 
' mentioned above plays an essentail role in the above unique phenomena, it also brings about some technical complexity 
of theoretical treatment, which hinders the elucidation of properties of weakly noncentrosymmetric superconductors. 
Moreover, the magnetic field induces the orbital depairing effect as well as the Zeeman effect mentioned above. To deal 
with these issues, we first analyze exactly effects of the Zeeman magnetic field on the quasiparticle energy, neglecting 
the orbital depairing effect, and, afterward, to take into account the orbital effect, we develope the quasiclassical 
method which extends the classical works by Eilenberger-^ and Larkin-Ovchinniko\*2i to the case with both the 
intra-band and inter-band pairs. A similar quasiclassical approach was considered before by Hayashi et al?^ Here, 
we obtain, for the first time, the explicit analytical solutions of the Eilenberger equations which encompass both 
the intra-band and inter-band pairs. Using them, we discuss behaviors of thermodynamic quantities under applied 
magnetic field, in which the above-mentioned unique features characterizing parity violated Cooper pairs appear. 

The organization of this paper is as follows. In Sec. II, we investigate the Zeeman field effect in weakly noncen- 
trosymmetric superconductor with the Dresselhaus type SO interaction, neglecting the orbital depairing effect, and 
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demonstrate that in the case of Ego A the Pauh depairing effect is anisotropic in the momentum space, leading the 
point-node-Hke anisotropic structure of excitation energy gap. This phenomenon should be important in the mixed 
state of type II superconductors. Thus, in Sec. Ill, we develope the quasiclassical method, taking into account the 
orbital depairing effect in addition to the above-mentioned anisotropic Pauli depairing effect. Based upon this method, 
in Sec. IV, we investigate the thermodynamic properties of the mixed state, and elucidate how the above features 
characterizing the parity- violated Cooper pairs appear in experimentally observable quantities. It is demonstrated 
that for a sufficiently small SO interaction, Eso < the specific heat as a function of magnetic fields exhibits a 
two-gap like behavior even when there is only a single SC gap. We shall also discuss the implication of our results 
for the recent experimental researches on the weakly noncentrosymmetric superconductor Y2C3i^'2£i^i2^ Summary 
is given in Sec.V. 



II. ANISOTROPIC PAULI DEPAIRING EFFECT 



The absence of the inversion symmetry is characterized by the asymmetric SO interaction, 

Ti-SO = " X! '^o(fe) • (Taa'cl^Cka', (1) 

where Cka (c^.g.) is the annihilation (creation) operator of an electron with momentum k and spin a. The components 
of <T = (CT^,cr^,cr^) are the Pauli matrices. Since we are concerned with the application to the cubic system Y2C3 
with the space group symmetry M3d, we assume the Dresselhaus interaction, Co{k) = (Cox, ^Oy, Cqz) — {kxiky — 

h.2\ jL (1.2 _ i.2\ jL (1.2 _ 1,2^)^) 31 

In the following, we consider the case of an s-wave state with the isotropic SC gap A, for which the features of 
weakly broken IS appear profoundly, as shown below. We ignore the admixture with triplet pairs, which is justified 
for Eso/E-p ^ 1. Our model Hamiltonian reads, 

^ = "^BCS + "^so , (2) 



Hbcs = ^£kci^Cka - ^[Acl^clf^^ + A*c_fcxcfc|]. (3) 

k,(T k 

In the following, the spherical Fermi surface, Ek = fc^/(2m) — Ep, is assumed for simplicity. In this section, we 
concentrate on the effects of the Zeeman interaction, expressed by — ■ (Tcya'C^^Cka-' with H = (0, 0, H^), 

leaving the analysis of the orbital depairing effect until the next sections. When the Zeeman term is added to the 
Hamiltonian the magnetic field induces the pairing between electrons on the SO split different bands, which 
gives rise to the Pauli depairing effect. The important observation is that in the case of £'so ^ A, Cooper pairs 
with k for which the spin degeneracy is not lifted by the SO interaction are more seriously affected by the Pauli 
depairing effect than electron pairs in the strongly SO split regions, leading to the anisotropic Pauli depairing in the 
momentum space. The anisotropic Zeeman effect in the strongly noncentrosymmetric case has already been discussed 
by several authors J^iii The unique point in the weakly noncentrosymmetric case is that this effect drastically changes 
the structure of low-energy excitations in the SC state, and yields the point-node-like structure of the single-particle 
excitation gap even in isotropic s-wave states. To demonstrate this, we calculate the single-particle excitation energy 
Ekf^ {fJ- = 1,2,3,4) by diagonalizing the above mean field Hamiltonian ([2]) with the Zeeman term expressed in the 




FIG. 1: (Color online) The structure of the single particle excitation energy gap on the spherical Fermi surface of the band Sk- 
for Hz = 0.0, 0.005, 0.01. 
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4x4 matrix form in the space spanned by the basis (c|.|, c_fe|, cj.|, C-ki)- The expHcit expressions for Ekfj. are given 
in Appendix A. In FIG.l, we depict the numericaUy calculated single-particle excitation energy gap i^gap for the band 
Ek- = £k — \c(Co{k) — fJ,BH\ in the case with A/Ep = 0.01, and a/Ep — 0.03. Here i?gap is defined by the magnitude 
of Ek3 for k on the Fermi surface satisfying Ek- = 0. It is seen that the point node-like structure develops, as the 
magnetic field increases. As a matter of fact, the excitation energy Ekf^ is not truely gapless, but the excitation energy 
gap £^gap is strongly anisotropic with the structure similar to the point nodes, reflecting the fc-dependence of the SO 
term even when the superconducting gap A is independent of k. For instance, the excitation energy gaps -Egap for 
kp 11 (001), kp II (111) and kp \\ (100) (equivalent to kp \\ (010)) are equal, and given by, 

- 4a"^ = EgT = ^J{mH.? + A^ - MBi?., (4) 

while for kp \\ (110), 

4ip"^ « \jyiH^) + - 2y^/(i/,)2 + ^,lH^^ (5) 

with f{H^) = 2a^k% + n^H^. It is easily checked that Egl°^ > E^l^K The positions of the gap minima coincide 
approximately with the zero points of Colk) at which the SO splitting vanishes. The quasi-particle excitations 
for kp II (001), (111), and (100) behave like Dirac fermions with mass gap given by eq.Q. It should be stressed 
that in the situation considered here the single-particle energy gap iS'gap does not coincide with the superconducting 
gap A, which is fc-independent. A similar anisotropic structure of the excitation gap also appears for the band 
ek+ = ek + \aCo{k) - hbH\. 

Although we use the Dresselhaus interaction in the present calculation, the point-node-like structure appears gen- 
erally for any forms of Co{k) which possess zero points on the Fermi surface. 



III. QUASICLASSICAL APPROACH FOR THE MIXED STATE IN THE CASE WITH BOTH 

INTER-BAND AND INTRA-BAND PAIRINGS 



The anisotroppic Pauli depairing effect discussed in the previous section is particularly important in type II super- 
conductors. Thus, in the following, we consider the orbital depairing effect as well as the Pauli depairing effect on the 
basis of the quasiclassical analysis. For this purpose, we consider the Green functions for quasiparticles defined on the 
SO split bands, G^t"'\x,x') = - {Trip ^{x)il}l{x')) , and F^^^^^x^x') = -{Tripl{x)i}l{x')) where ip^,{x), ip^^x) are the 
field operators for quasiparticles in the /i-band corresponding to the energy in the normal state = Efc -I- ^.a\C^{k)\ 
with /i = As mentioned above, the inter-band Green functions plays important roles in addition to the intra-band 
Green functions. Fourier transforming G('^'^) and F'^^^'^\ we introduce the quasiclassical Green functions defined by 



C/(fc,r,e„) 



/( + + ) g( + + ) /(- + )t 
g(- + ) ^f{- + ) g( — ) 
V/( + -)t -gi + -) /(-)t 



( + -) _/( + -) \ 



(- + ) 

-ni — ) 



(6) 



with5('^-)(fc,r,£„) = / i|i.G('''')(fc,r,e„), Z^^'') (^, r, £„) = / (fc, r, £„), and gf^'') = g(^'')(-fc, r, -e„). Here 

k is the momentum conjugate to the relative coordinate x — a;', r is the center of mass coordinate. A; is a unit vector 
parametrizing the direction of momentum fc, and £„ is the fermionic Matsubara frequency. Hereafter, matrices A 
represent 2x2 matrices defined in the space spanned by the basis (ipf^jipl) where ^ = i' = ± or fi = —v = ±, 
and matrices B represent those in the two-dimensional space spanned by the band indices — . Using the standard 
method^^i^iiH we find that Q satisfies the Eilenberger equation in the clean limit. 



iv ■ —g + [wr" - M A, = 0, 
or 



(7) 



where u = 2i£„ + v ■ ^A with A a vector potential, and t^' — 1® . The 4x4 matrix M is defined by M{k) = 

fiBH\. The 



CT^ {aL{k)), where L{k) = l(|£(fc)| + |£(-fc)|)/2 + (T"(|£(fc)| - |£(-fc)|)/2, and a|£(fc)| = \aCo{k) 
matrix gap function A(fc,r) is 



A(fc,r) = 



A+(fc,r) A2(fc,r) 
-A2(fc,r) A„(fc,r) 



(8) 
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FIG. 2: (Color online) Specific heat at T = O.lTc divided by the normal state value Cat as a function of the normalized magnetic 
field h = ^B-ffz/Ao for q/Aq = 1.0 (solid line), a/Ao = 0.0 (dotted hne), and a/Ao = 10.0 (dashed line). 



with A±(fc,r) = i(Tj,ReA±(fc,r) + io-3;ImA±(fc, r), A2(fc,r) = aa;ReA2(fc, r) - o-j^ImA2(fc, r), A±(fc,r) = A(r)s±(fc), 
A2(fc,r) = A(r)s2(fc), and, 

s±Ck)=-[Mk)i-i-k)+^-m+i-k)]vTik) (9) 



,S2(fc) - e+(fc)e+(-fc) - (10) 



where ^±{k) = y {1± Cz{k)\C{k)\)/2, and ?7±(fc) = —{Cx ± iCy)/ ^ + Here, we omit the corrections to the 

Fermi velocity from the SO interaction, which are not important for Eso/E-p <C 1. In the case with strong SO 
interaction a 3> A, the contributions from the inter-band Green functions to eq.([7]) are negligible, and the two bands 
are decoupled. This simplified case was previously studied by several authors i^i^^^ 

In spite of the complications raised by the existence of both the inter-band and intra-band Green functions, the 
analytical soutions of ([7]) based on the Pesch type approximation are possible. In this approach, we assume the 
Abrikosov lattice solution for A(r) and the uniform magnetic field in the system, and replace the normal Green 
functions g''^'^\ ^(a"^) with the spatial averages over a unit cell of the vortex lattice, {g'^^"^), (g*-'"^''), retaining only 
the spatial variation of and Utilizing the normalization condition Q ■ Q = 1 and the relations tr[Q] — 0, 

g{A )* _ g{ h)^ which are derived from we find that the quasiclassical Green hmctions are given by the solutions 

of the coupled algebraic equations, 

(T.I/— ± 

-vos2{{g^""^^^)+v{g^^---^))\'^-l, (A* = ±), (11) 



{g^+-^) J2 '^[Yis<rW'"''>) ~ Yieir "W"''')], (12) 

where sit""^ = e„ + - u\C'\), C = £(-fc), P{e) = ^dYie)/ds Yjin) = V^A^UnW{2iUnSn) with A^ = 

(A^(r)), W{z) = e"^ erfc(— iz), and w„ = A sgne„/usin0 with A = hc/2eHz and 6 the polar angle of fc, and 
C = 2is2S+a\C'\/D with 

D = ia'\C\\C'\ + 2^\s+\M\C\ + \C'\)Y,^^Yielr^) 

-2^sla{\C\ - \C'\) E f^Y{elr^). (13) 

Since the explicit expressions of the solutions are lengthy, we will present them in Appendix B, and, instead, concen- 
trate on the discussion on the results obtained from them in the following section. 
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FIG. 3; (Color online) Log- log plot of the density of states versus energy for a/Ao = 1.0. The density of states for a = 0.0 and 
ft = 0.11 is also shown for comparison (broken). The dotted line is . 

IV. SPECIFC HEAT COEFFICIENT AND DENSITY OF STATES 

Using the analytical solutions for we calculate the specific heat and the density of states of quasiparticles. For this 
purpose, we determine the field-dependence and the temperature-dependence of the spatially averaged gap function 
A{Hz,T) = -y/ (A2(r)), by solving the quasiclassical BCS gap equation, 

(A2(r)) = AoT^^(A(r)/Ti(fc,r,£„)), (14) 
" k 

fn = + - - i-^+f^—^ (15) 

where = S,ai—k). Ao is a dimensionless coupling constant. To compare our calculated results with the experimental 
observations for Y2C3 later, we tune the values of Aq and the cutoff for the frequency sum Sc so as to realize Tc — 18 
K. (e.g. Aq — 0.2703, Sc — 700 K.) In this calculation, an important parameter is ao = \/ hb^c/ (2eAo)/(7r^o), which 
is an inverse of the coherence length ^0 normalized to be dimensionless. Here Ao is the gap function for = 0. For 
Y2C3, ao « 0.234. The specific heat in the SC state is 

Cs= I de Ds{e), (16) 

J 4T2cosh2^ ^ ' ^ ' 

where the density of states for quasiparticles Ds{e) is given by Ds{e) = Dn{0) J^k X]^=± Iin(5''^^''(A;, r,e + iS)), with 
D]^{0) the density of states in the normal state. In FIGlH we present the calculated results of the specific heat as 
a function of the normalized magnetic field h = fiBHz/^o for the temperature T/Tc — 0.1 and several sets of the 
parameters a/Ao and ao- It is found that in contrast to the case with strongly broken IS where the SC state is 
quite robust against the Pauli depairing effect, in weakly noncentrosymmetric systems more low-energy excitations 
are induced by the magnetic field than in the case with inversion symmetry a = 0. We interpret the origin of this 
behavior as the existence of the field-induced nodal excitations mentioned above. The emergence of the point-node- 
like excitations is more clearly observed in the energy-dependence of the density of states for quasiparticles, which 
is plotted in FIGl3]for a/Ao = l-O- As the magnetic field increases, the density of states exhibits the power-law 
behavior Ds{e) cx e^. 

When the SO coupling a is sufficiently smaller than Aq, another remarkable feature appears in the field-dependence 
of the specific heat. As shown in FIGIH for the low magnetic fields h « 0.1 ~ 0.2, a shoulder-like structure of Cs{Hz) 
similar to a two-gap behavior appears, though there is only a single SC gap A(r). The origin of this behavior is 
understood as follows. For a < Aq, there are two different types of the Pauli depairing effect; i.e. (i) one due to the 
generation of the inter-band pairing correlation which, instead, suppresses the intra-band pairing, and (ii) the other 
caused by the asymmetric deformation of the Fermi surface. The former is analogous to the usual Pauli depairing effect 
in centrosymmetric superconductors. The latter effect is inherent in noncentrosymmetric systems. The crucial point 
is that although the former exists for any finite magnetic fields, the latter is effective only for small fields ^jlbHz < a, 
and is suppressed for larger fields finHz > a. As a result, the character of the Pauli depairing effect changes around 
Hz ^ a/ in the case of a < Aq, yielding the shoulder structure of the specific heat coefficient as demonstrated in 
FIGIH We would like to stress that this effect is caused by the momentum dependent spin orientatin of Cooper pairs 
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FIG. 4; (Color online) Specific heat as a function of the normalized magnetic field h for small o/Aq. oo = 0.234. T = O.lTc 

characterizmg parity violation. Three remarks are in order. (1) For sufficiently small a, e.g. a ~ 0.05, Cg exhibits 
a hump rather than a shoulder, in marked contrast with conventional two-gap behaviors. (2) As a increases, the 
position of the shoulder shifts to larger h regions, though its structure becomes obscure since the orbital depairing 
effect dominates for high magnetic fields. (3) It should be cautioned that the Pesch approximation is not valid for 
magnetic fields much lower than _ffc2-— However, the shoulder structure of Cs at h ^ stems from the Zeeman 
effect rather than orbital depairing effects. Thus, the above results are applicable to systems with a sufficiently large 
value of the Ginzburg-Landau parameter k, where small magnetic fields can penetrate deeply into the SC regions. 

We, now, discuss the implication of the above results for the experimental observations of the weakly noncentrosym- 
metric superconductor Y2C3. This system is almost in the London limit with n > 10.^'^'?- According to the recent 
LDA calculation, the averaged magnitude of the SO band splitting is roughly ~ 0.01 eV,'^2, which is of the same order 
as the SC gap ~ 30 Kf22 Thus, our analysis is applicable. A remarkable experimental observation for Y2C3 is that the 
field-dependence of the specific heat exhibits a small shoulder structure for ^ 8 ~ H^il'^ T at T = 2.6 Ki^ Al- 
though this behavior was interpreted as the indication of the existence of two SC gaps with different magnitudes 
it can be also, alternatively, explained by assuming the realization of the unique effect associated with parity violation 
as demonstrated in FICU) A possible test for our scenario is to investigate the field dependence of the nuclear spin 
relaxation ratCj^- which should exhibit a gap energy scale different from that observed in the specific heat. 

V. SUMMARY 

We have shown that in weakly noncentrosymmetric superconductors the Pauli depairing effect is anisotropic in the 
momentum space, inducing the point-node-like structure of the quasi-particle excitation energy gap even in isotropic 
s-wave states. This effect is caused by the competition between the asymmetric SO interaction and the Zeeman 
magnetic field in the superconducting state, and yields unique low temperature behaviors of thermodynamic quantities 
quite different from those of conventional s-wave superconductors. Also, by using the quasiclassical method, we have 
demonstrated that the magnetic field dependence of the specific heat exhibits a multi-gap-like structure for sufficiently 
small SO interaction even when there is only a single gap. These effects are associated with the momentum-dependent 
spin orientation of Cooper pairs which characterizes parity violation. Thus, our results reveal the unique aspects of 
parity violated Cooper pairs inherent in weakly noncentrosymmetric systems. We have also discussed that our findings 
may be relevant to the recent experimental observations for Y2C3. 
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APPENDIX A: QUASIPARTICLE ENERGY 



The single-electron energies of the Hamiltonian ([2]) with the Zeeman term — ^J.BHzU^„ic\^Cka' are given by 

the solutions of the following eigen value equation, 



Efc — MB^^z + -Coz — X £02; — i^Qy 

-Sk + + Az - a; -A* 

Cqx + i-Coy -A Efe + /zb^^z - -Coz - a:^ 

A* 



A 




0. (Al) 



Here the matrix is defined in the space spanned by the basis (c^p c_fe|, c^^, c-fcj^). The explicit solutions of (|A1|) are 

/ 



V 



-2Y 



1 1 



<5| +(5i - |X 



114 



(A2) 



£^/c2 — - 



1 1 



& ' - -X - 
- 3 



2Y 



al + a: 



-1^ 



3 



(A3) 



-Efcs — - 




(A4) 



EkA — - I — 



1 1 



ai + &i- -X 



\ 



2Y 



al — at X , 



(A5) 



where a 



^ ^ e^'a'^Q-l., and 



with 



1 . 1 



a± = ^{j;^{2X^ + 36XZ)-AXZ + 



± ^/(^(2X3 + 36XZ) AXZ + y2)2 _ _l(i2Z + A2)3) 



(A6) 



X = -2{el + a^\C^{k)\^ + /^fff + A^), 



(A7) 



Y = 8pLBHzaCoz£k, 



(A8) 



Z = (4-a^|/:o(A:)|^-MB^fz)' + A* + 2AVfe-4Ai^F>^£^,. 



(A9) 



Since there is a term linear in x in eg. ljAip with the coefficient given by (jASp for Hz 7^ 0, the particle-hole symmetry 
is broken by the Zeeman magnetic fieldj ^^'^^ 



APPENDIX B: SOLUTIONS FOR QUASICLASSICAL GREEN FUNCTIONS 

In this appendix, we present the explicit solutions for the spatially averaged quasiclassical Green functions sat- 
isfying egs. ljlip and p2p . which are derived from the Eilenberger equations ([7]) combined with the normaliza- 
tion condition Q ■ Q — 1. Using the decoupling approximation for the spatial average, (g^^*^^ ^) « {g^'^'^^)^, 
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(g(/''^)/(«A)^^ ~ {gi^'•')){f'^^)A), etc., we obtain the expressions for (g^'''')), 

{g(++\p)) = -isgne„[l-P(4++))|s+(fc)|2{l-r(p) (Bl) 

- 2a{\C'\~\C\)s,{k){a{p) + b{p)r{p))f 

+ 4a'\C'nn{k)\\a{p) + b{p)r{p)f 

- P(4-)){s2(fc)(l + r2(-p)) 

+ 2a{\C'\ + \C\)\s+ik)\'{a{p) + 6(p)r(p))}^]-^ 



where 



{9^--\p)) = -rip)r2{pW++\p)), (B2) 



a{p) = '-^[Y{1) - Y{3) + M-p){Y{2) - Y{3))], (B3) 

Kp) = -'-^[Y{1) - Y{A) + r2{-p)iY{2) - y(4))], (B4) 

Mp) = -^[A_+8a3(|£'|_|£|)(|£|2 (B5) 
- \Cf)\s+{krs,{k)4Y{l)], 

A± = Y{1) - 2a{\C'\ - \C\)s2{k)coY{l) (B6) 

- 4a|i:|s2(fc)cor(3) 

+ 4a'{\C'\ + \C\)\s+{k)\'cl{Y{4)\C'\ - Y{3)\C\) 

± 8a2(|£|2 _ |£f )co|s+(fc)|2 

X [c±{F(l) - 2a|£|y(3)coS2(fc) 

- a{\C'\ - \C\)s2{k)coY{l)} 

± acls2{k){Yi3)\C\TY{4)\C'\)] 

c+ = ^[F(2) - + y(4) - y(3)] (B7) 

c_ = !f^[r(4) + r(3)-2r(2)] (bs) 

co = if^[r(l)+y(2)-2y(3)] (B9) 

r-(p) = ^[-^ + sgn eWS^-4RQ], (BIO) 

g = y(l) - 2a{\C'\ - \C\)s2{k)Y{l)a{p) (Bll) 

- 2a|£'|y(4)(l + r2(p))s2(fc)a(p) + a{pfh{p), 

R = -Y{l)-2a{\C'\-\C\)s2Ck)Y{l)b{p) (B12) 
+ 2a|£'|r(4)(l + r2{p))s2{k)b{p) + b{pfh{p), 
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S = -2ai\C'\-\C\)s2ik)Y{l){a{p)+bip)) (B13) 
+ 2a|£'|r(4)(l + r2(p)).S2(fc)a(p) 
- 2a|£|r(3)(l + r2{-p))s2ik)b{p) + 2a{p)b{p)h{p), 



with 



h{p) = 4a'{\C'\ + \C\)\s+{k)\' (B14) 
X (y(4)|£'|-y(3)|£|). 

Herep = (fc,e„), = Yiei+^^), Y{2) = Y{3) = Y{el'~^), and Y{4) = Y{ei+-^). 

Substituting the above expressions for (g*-''^') into cq. (I12|) . we obtain the off-diagonal components of the normal 
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{g^-^\^p))--^^y-^\p))- (BI8) 
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